Quantization from Hamilton-Jacobi theory with a random constraint 
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We propose a method of quantization based on Hamilton-Jacobi theory in the presence of a 
random constraint due to the fluctuations of a set of hidden random variables. Given a Lagrangian, 
it reproduces the results of canonical quantization yet with a unique ordering of operators if the 
Lagrange multiplier that arises in the dynamical system with constraint can only take binary values 
±h/2 with equal probability. 
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I. MOTIVATION 

In the previous work we have developed a new 
method of quantization of system of spin-less particles 
based on a specific modification of the classical dynam- 
ics of ensemble of trajectories parameterized by an unbi- 
ased (hidden) random variable. The fluctuations of the 
random variable is then characterized by the Planck con- 
stant. Given a wide class of classical Hamiltonians, the 
quantization is done by assuming rules of replacement of 
(deterministic) c-number by (stochastically parameter- 
ized) c-number to be applied to the classical Hamilton- 
Jacobi and continuity equations generated by the clas- 
sical Hamiltonian. Assuming a specific type of distri- 
bution of the random variable, the modified Hamilton- 
Jacobi and continuity equations are then written into the 
Schrodinger equation with a unique quantum Hamilto- 
nian, thus is free from the operator ordering ambiguity 
of the canonical quantization [2j. Moreover, one can al- 
ways identify an effective velocity field which turns out 
to be equal to the actual velocity field of the particles in 
pilot- wave theory [3j . 

In the present paper, we shall show that the rules of 
replacement heuristically postulated in Ref. 1] can be 
derived from the Hamilton-Jacobi theory with a random 
constraint. The fluctuations of the constraint is assumed 
due to the presence of some background fields whose de- 
tail interaction with the particles are not known. The 
existence of background fields is also assumed in the Nel- 
son stochastic mechanics to explain the origin of quantum 
fluctuations as time symmetric Brownian motion [4j . We 
shall then show that the hidden random variable postu- 
lated in Ref. [l[ can be identified as the Lagrange multi- 
plier that arises in the corresponding dynamical system 
with random constraint. 



II. HAMILTON-JACOBI THEORY WITH 
RANDOM CONSTRAINT 

Let us first review the Hamilton-Jacobi theory of clas- 
sical dynamics. For simplicity, below we shall discuss 
the case of system with finite degrees of freedom. Let 
us denote the corresponding Lagrangian of the system as 



L(q, q) where q = {qi}, i = 1,2,..., runs for all degrees 
of freedom, is the configuration coordinate, t is time and 
q = dq/dt = {qi} is the velocity. The canonical conjugate 
momentum is defined as 



~ l dqi 



(1) 



If the Lagrangian is not singular, det(9 2 L/dqidqj) ^ 0, 
which is assumed to avoid unnecessary complication, 
then the above equation can be solved in term of q as 



(2) 



Now let us define action as I = J L(q, q)dt. The actual 
trajectory connecting two spacetime points is then ob- 
tained by extremizing 7 with respect to variations of (q, q) 

for a pair of fixed ends SI_ — S f^^-j L.(q, q)dt — 0. This 
is the Hamilton's principle of stationary action which 
leads directly to the Euler-Lagrange and Hamilton equa- 
tions. 

Another way to solve the above variational problem is 
through the Hamilton-Jacobi method [5]. First we con- 
struct a one-parameter family of hypersurfaces S_(q; t) = 
t, where r is some parameter labeling the hypersurfaces, 
so that any point in the configuration space belongs only 
to one of the surfaces. Moreover a trajectory crosses each 
hypersurface only once and is nowhere tangent to it. r 
is thus a function of time, r = r (t) , so that one has 



. . dr dS _ . „ „ 



(3) 



By construction we have A / 0. The Hamilton-Jacobi 
method to get the actual trajectory then proceeds as 
follows. Let us consider two hypersurfaces S_ = r and 
S_ = r + dr, separated by infinitesimal dr. Then given 
any point on the hypersurface S_ = r, the actual tra- 
jectory that passes through this point and reaches the 
surface of 5 = r + dr is the one with a velocity q that 
minimizes dj_/dr Hence, the actual velocity must 
solve the following necessary Hamilton-Jacobi condition: 



d_ 

dq 



7, \drJ 



(4) 
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with fixed dr. Since dl/dr = L/A and A ^ 0, Eq. © 
then reduces into 



(5) 



where we have used dA/dtji = d qi S_ from Eq. (J3]). 

Next, assuming that L/A — 8(t), where 9 is a function 
only of t, then one can reparametrize S_ so that L = A 
In this case, one gets 



dL 



p -^wr°^- (6) 

Inserting this into Eq. one thus has 



qi(q,d q S;t). 



(7) 



Finally, substituting Eq. (JT)) back into Eq. ([3]), recalling 
L = A, one obtains the Hamilton- Jacobi equation 

d t S + q(q, q S; t) ■ q S - L(q, q(q, q S; t)) = 0. (8) 

One can further show that the Hamilton- Jacobi equa- 
tion of ([5]) is equivalent to the Euler-Lagrange and the 
Hamilton equation [5] . In contrast to the latter two equa- 
tions, the Hamilton- Jacobi equation describes a congru- 
ence of trajectories whose velocity field is given by Eq. 
([7]). A single trajectory is then obtained if one also fixes 
the initial configuration of the system. Another impor- 
tant feature of Hamilton- Jacobi theory is that it imposes 
a local condition to the actual trajectory, that of Eq. ((4]). 
By contrast, Hamilton's principle gives a global condition 
to the actual trajectory. 

Now let us assume that the system under consid- 
eration depends on a set of hidden random variables 
£ = (£i , £2 j ■ • ■ ) whose dynamical origin is not known. 
This for example might be due to the presence of back- 
ground fields whose detail interaction with the particles is 
not known resulting in a stochastic motion of the latter. 
Single event is thus inherently random. Hence, one can 
only make prediction concerning an ensemble of copies of 
the system. Let us denote the joint-probability density 
of the fluctuations of q and £ as ft(g, £; t). The marginal 
probability densities of the configuration of the system q 
and £ are then given as 

p(q;t) = J d&(q,Z;t) & P(0 = J dqSl{q,t;t). (9) 

Here we have assumed that the probability density of £ is 
independent of time. Let us further assume that there is 
a random velocity field v(q, £; t), of the same dimension as 
the system, so that ft(g,£;i) has to satisfy the following 
differential equation: 



G(q,q]tt,v) 



d(lnft) 
~dt 



-dq-V = 



ft 



d q n 



no; 

The above constraint might be interpreted that the ve- 
locity divergence gives the only source of (local in config- 
uration space) change of entropy. 



Next let us construct a new family of hypersurfaces 
S(q,£,it) — T so that for a fixed value of £, any point in 
configuration space belongs to only one of the hypersur- 
faces. Yet, a single point can belong to more than one 
hypersurfaces with different values of £. Namely S(q, £; t) 
is now fluctuating due to the fluctuations of £. Further, 
let us assume that for a fixed value of £, a trajectory 
satisfying Eq. (|10|) crosses a hypersurface only once and 
nowhere tangent to it. r is thus a function of t and £, 
t = r{t 1 £). Now, let us again consider two hypersurfaces 
S = t and S = t + dr separated by infinitesimal dr. Let 
us assume that within this interval, £ is fixed. Then, as in 
the case with no constraint discussed before, let us pos- 
tulate that given any point on the hypersurface S = r, 
the actual trajectory that passes through this point and 
reaches the hypersurface S = r + dr is the one with a 
velocity q that minimizes dj_/dr and satisfying the con- 
straint of Eq. (fTU|) . We thus have to solve Eq. ((U) with 
the constraint of Eq. (flUl) . Using the Lagrange method, 
this problem implies the following necessary condition: 



dqi \drJ dqi 



0. 



(11) 



where A = A(£) is the Lagrange multiplier. 

Notice that since the constraint is fluctuating ran- 
domly due to the fluctuations of £, then the Lagrange 
multiplier also depends on the value of the hidden vari- 
ables thus is inherently random. Let us denote the prob- 
ability density of A as -P(A). As is clear from Eq. (jlll) . 
£ appears explicitly in the equation only through A(£). 
Accordingly, we shall regard A as the effective hidden 
random variable and use it in place of £. Hence, for ex- 
ample, we shall write 51(g, A; t) instead of Vl{q,^\ t) and so 
on. Further, by construction £ thus A in general depends 
on the configuration space and time. Later we shall as- 
sume that the derivatives of A with respect to space and 
time are negligible as compared to that of S. The value 
of A can be obtained by inserting Eq. (|11[) back into Eq. 
(|TU| . A thus depends on and v. Below we shall go 
the other way around. Namely, we shall assume A with 
a specific statistical properties and look for a class of 51 
and v that satisfy Eqs. (|T0|) and (fTT|) . 

One can then proceed as before to arrive at the follow- 
ing pair of equations [f| : 



dL . Q . 

dqi dq 



0. 



(12) 



where L = L(q, q, A; f2, v) is an extended Lagrangian de- 
fined as 



L(q, q, A; ft, v) = L{q, q) - \G{q, q; ft, v) 

= L(q,q ] t)-X(^ + ^-q + d q -v). (13) 

Equation (TT2"]) has to be solved together with the con- 
straint of Eq. PH|) . From the left equation of (TT21 . using 
Eq. (flUl) . one gets 

dL 8 ft 

Pi(q, q, A; t, ft) = — = p. (q, q; t) - A^- = d m S. (14) 
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Again, assuming that the unconstrained Lagrangian L is 
not singular, the above equation can be solved in term of 
q to give 



q(q, d q S + A^;t). 



(15) 



This velocity field implies the following continuity equa- 
tion conserving the probability: 



dq ■ (q{q, d q S + A 



d q n 



o. 



(16) 



Further, inserting Eq. (|15p into the right equation of 
(|12p. one obtains the following (modified) Hamilton- 
Jacobi equation: 

-L(g, <j(g, q S + A^; t), A; Q, w) = 0, (17) 

where we have used the left equation of (fT2")l . 

We have thus Eqs. (TTB]) and (|17l) that have to be solved 
in term of fi and S subjected to the condition of Eq. (pT0|) . 
They are coupled to each other since q now depends on S 
and O. We shall show in the next section by taking a con- 
crete example that Eqs. (IT6|) and (jXTJ) can be combined 
together to give a single partial differential equation for 
S and Q Q- The combined equation has to be solved 
with the condition of Eq. (fTU|). To do this, we thus have 
to express v that appears in both equations as function of 
S and fl. On the other hand, to be meaningful, the con- 
straint has to be consistent with the dynamics. Namely, 
if initially ft satisfies the constraint of Eq. (fTU|) then it 
must be so for any time as fl is evolved by the dynamics 
through Eq. (|16j) . This naturally implies that the ran- 
dom velocity field of the constraint v has to be related to 
velocity field of the dynamics q. 

Now let us assume that v(q, A; t) is given as 



Vi(q,X;t) 



qi(q,X;t) +q l {q,-\;t) 



Vi(q,-X;t). (18) 



v(q, A; t) is thus uniquely determined by the choice of 
unconstrained Lagrangian L. Moreover, let us further 
assume a class of solutions satisfying the following sym- 
metry relations: 



S(q,X;t)=S(q,-X;t)+S (X) & tt(q,X;t) = Q(q, 



-A;t), 

where So (A) is independent of q and t. The former can be 
done by choosing an appropriate parameterization of the 
hypersurfaces, namely r(A,t) = r(— A, t) + So (A). More- 
over, the latter implies that A is an unbiased random 
variable 



P(A) = J dqQ(q,X;t)=P(-X). 



(20) 



Let us show that the above choice of random velocity 
field v generates a constraint that is consistent with the 



dynamics. To see this, first, notice that taking the case 
when A is positive in the constraint of Eq. (ITUl) add to it 
the case when A is negative and divided by two, imposing 
Eqs. (JTSJ) and (fI5]l. one gets 



dtQ. + d q ■ (vQ) = 0, 



(21) 



which is just a continuity equation. Hence, in this case, 
the constraint is just a probability conservation equation 
generated by random velocity field v. On the other hand, 
from Eq. (fT6|) , taking the case when A is positive add to 
it the case when A is negative and divided by two one 
gets, by virtue of Eqs. ([I5J and (TK)1) . 



d t n + d q ■ (vft) = 0, 



(22) 



which is the same as Eq. (|2"Tj) . Hence the constraint is 
consistent with the dynamics, as expected. 



III. "EFFECTIVE" PILOT- WAVE MODEL 

For simplicity, let us apply the general formalism de- 
veloped in the previous section to ensemble of system of 
single particle of mass m subjected to external potentials. 
The unconstrained (classical) Lagrangian then takes the 
form 

L(q,q) = ±mq 2 +A(q) -q-V(q), (23) 



so that one has 

p. = mqi + Ai. 
The classical Hamiltonian reads 

2m — 

Inserting Eq. into Eq. (fbfj) one gets 

d Qi S Ai A d„.n 

Qi = — + 

m m m 11 

Hence, Eq. (TTB1) becomes 



(24) 



(25) 



(26) 



1 



A 



d t n + — d v ■ Ud„s - A)n) + — o 2 „n = o. (27) 

m m q 

Next, substituting Eq. PS| into Eq. (|18p and imposing 
the assumption of Eq. (fl"9|) , v is given by 



(28) 



Vi(q,X;t) = —(d qi S-Ai). 



Inserting Eqs. (p ]l , ([23 ]l . ([26 )l and ([28 ]l into Eq. $T7$ one 
obtains 

_ {d q S-Af 2X 2 d 2 n R 

Oto + 



2m 



V 



m R 



+ Ti ( dtil + ^ ' ( {d « S ~ A)n ^ + ^) = °' (29) 
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where we have defined R = V" and used the following 
identity: 



ld qi Vld qj Q, _ ld qt d qj il d qi d qj R 
4 ^ ~~ 2 U R ' 

Substituting Eq. ([2"7]l. Eq. (gHJ) then reduces into 



(30) 



^ -T ^ 0. (31) 



2to 



Moreover, in this case, inserting Eq. (f2"5]l into Eq. ([2~T]) . 
the constraint then reads 



+ • ((9,5 - A)0\ = 0. 



(32) 



We have thus pair of coupled Eqs. (|3Tj) and (|32l) pa- 
rameterized by the fluctuating Lagrange multiplier A(£). 
Now, since A is non-vanishing, one can define the follow- 
ing complex- valued function: 



*(g,A;t) = J Rexp(^j5). 



(33) 



It differs from the Madelung transformation in that S 
is divided by 2|A| instead of H so that one has p(q; t) — 
J d\n = J dA|*| 2 . The pair of Eqs. ^ and ^ can 
then be recast into the following modified Schrodinger 
equation parameterized by the fluctuating Lagrange mul- 
tiplier: 

i2\\\d t V = l-(-i2\X\d q - A)H + VV, (34) 

where we have imposed the assumption that the space 
and time derivatives of A are negligible as compared to 
that of S. 

Let us proceed to assume that 51 is factorisable as 



Ct(q,X;t) = p(q,\X\;t)P(X), 



(35) 



with J dqp{q, |A|; t) — 1 for arbitrary value of A. This 
guarantees that O is correctly normalized J dXdqil = 1. 
Moreover, let us assume that the Lagrange multiplier can 
only take binary values A(£) = ±h/2 with equal proba- 
bility 

P(A) = ~<5(A- ft/2) + 1,5(A + ft/2). (36) 



In this case, Eq. (j3"4"]l becomes 



(37) 



where the wave function ^Q(q;t) is given by 

*q(95 <) = Vpfe V2; *) exp (|%, ±ft/2; t)). (38) 

Hence p(q; t) — \^fg(q; t)\ 2 holds by construction, and the 
phase is given by Sq{q;t) = S(q,±H/2;t). 



The quantum mechanical Schrodinger equation is thus 
reproduced as a specific case of the present statisti- 
cal model when the Lagrange multiplier A is an unbi- 
ased binary random variable which can only take values 
A(£) = ±ft/2. To this end, it is interesting to mention 
Ref. Q which showed that the master equation of a par- 
ticle moving with a fixed velocity, imposed to a random 
complete reverse of direction following a Poisson distri- 
bution, can be written into Dirac equation (in the same 
way that the Schrodinger equation is connected to the 
dynamics of Brownian motion) through analytic contin- 
uation. Note also that while A is a binary random vari- 
able, it is a function of the true hidden random variables 
£ which in turn may t ake continuo us values. For example 
one may have A = Vsi+Cl+W = ^ft so that £ lies on 
the surface of a ball of radius ft If we divide the surface 
of the ball into two, attribute each devision ±ft, respec- 
tively, and assume that £ moves sufficiently chaotic, then 
one will have A = ±ft with equal probability. 

Moreover, in this case, the velocity field v that gener- 
ates the constraint then becomes 

Vi(q,±h/2;t) = —(d qi S(q, ±ft/2; t) - Ai) 
m 

= ^(d qi S Q (q;t)-A i ) = v i (q;t). (39) 
Keeping this in mind, Eqs. ([2"Tj) or (|3"2")l then reads 



d tP +^d q -({d q S Q -A)p) =0, 



(40) 



where we have used Eq. (|55|) . It turns out that v defined 
above is numerically equal to the actual velocity of the 
particle in pilot- wave theory [3]. This is also equal to 
the "naively observable velocity field" reported in Ref. 
Q , obtained using the notion of weak measurement [Tpj 
within the standard interpretation of quantum mechan- 
ics. 

We have thus an "effectively" similar picture with 
pilot-wave theory in the sense that the particle always 
possesses definite position and momentum and further it 
moves "as if" it is guided by the wave function so that 
the "effective velocity" v is given by Eq. (I3"9l . Hence, we 
can conclude that the statistical model developed in the 
present paper will reproduce the statistical wave-like in- 
terference pattern in slits experiment and tunneling over 
potential barrier It is then tempting to further in- 



vestigate in the future whether the model can lead to a 
description of quantum measurement which solves the in- 
famous measurement problem. Note that to discuss the 
problem of measurement one needs to consider the time- 
irreversible process of registration which involves realistic 
description of apparatus and bath with large degrees of 
freedom [lS^ . 

However, in contrast to pilot-wave theory, as is clearly 
shown in Eq. (I31)l) . v is not the actual velocity of the 
particle, but the average of two actual velocities corre- 
sponding to Lagrange multiplier equal to h/2 and —ft/2. 
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Hence, while pilot- wave theory is strictly determinis- 
tic, the present model is inherently stochastic. Further, 
in this statistical model, rather than being postulated, 
the Schrodinger equation and a unique guidance relation 
emerge naturally by imposing a random constraint to 
the Hamilton- Jacobi condition. It is also evident that 
in contrast to pilot-wave theory, the wave function in 
the model is not physically real. It is just an artificial 
mathematical tool to describe the dynamics and statis- 
tics of the ensemble of trajectories, and by construction 
the Born's statistics p(q;t) = \^Q(q;t)\ 2 is valid for all 
time. The so-called quantum potential given by the last 
term of Eq. (|3lT) with A = ±7i/2, which is argued by 
pilot-wave theory to be responsible for all peculiar quan- 
tum phenomena Q, is generated by the local change in 
configuration space of the entropy part of the constraint. 
Note however that since the Schrodinger equation is time- 
reversal invariant then the total change of entropy should 
be vanishing. Finally while pilot- wave theory can deal 
with a single trajectory (since the wave function is as- 
sumed to be physically real satisfying the deterministic 
Schrodinger equation), the present model strictly con- 
cerns the dynamics of an ensemble of trajectories as in 
Nelson stochastic mechanics. 



IV. QUANTIZATION WITH UNIQUE 
ORDERING AND DIRECT PHYSICAL 
INTERPRETATION 



where H_(q,p) = q ■ p — L has the same form as classical 
Hamiltonian. Next, applying the left equation of (|41|) . 
the velocity field of Eq. ([T5j) also becomes 



q{q,p;t)\ P =d q s 



dH 
dp 



P=d q S 



(44) 



so that the continuity equation of Eq. (fT6|) now reads 

dH 



dp 



P=d q S 



0. 



(45) 



Notice then that Eqs. (1431) and (|4"5")) take the same 
form as the Hamilton- Jacobi equation and the continuity 
equation of classical mechanics given a classical Hamil- 
tonian H_(q,p). Hence, given the classical Hamiltonian, 
to get the Hamilton- Jacobi equation of (|17l) and con- 
tinuity equation of (|16[) based on which we derive the 
Schrodinger equation, we can first develop the corre- 
sponding classical equations of (|4Uf and (1431) . and apply 
the rules of Eq. (|41[) . The remaining task is to express v 
defined in Eq. ([18)) in term of the classical Hamiltonian 
H . To do this, one can see that if the classical Hamil- 
tonian is at-most-quadratic in classical momentum, then 
q is a linear function of p = d q S_ = d q S + X(d q Q/Q). In 
this case, evaluating Eq. ([T8"l) and taking into account 
Eq. (TH)1) one then obtains 

q(q, 8 q S + A^; t) + q(q, 8 q S - A^; t) 



We have shown in the previous section by taking an 
example of particle in external potentials how to develop 
from a given classical Lagrangian a Schrodinger equa- 
tion with unique quantum Hamiltonian. It can thus be 
regarded as to provide a method of quantization of classi- 
cal system given its classical Hamiltonian. In contrast to 
the canonical quantization which is formal-mathematical, 
the method presented in the previous section is based on 
Hamilton- Jacobi theory in the presence of random con- 
straint so that the quantum-classical correspondence is 
physically kept transparent. 

To develop further formal comparison with canonical 
quantization, now let us extract some rules which can 
be applied directly given a classical Hamiltonian. First, 
let us define a scalar function S_ so that dS_ — dtSdt + 
d q Sdq = Ldt. Then, recalling that dS — dtSdt+d q Sdq = 
Ldt, from Eq. (IT3"1) . we have the following pair of rela- 
tions: 

d q S = 8 q S + A^, d t S = d t S + A^ + \d q ■ v. (41) 

On the other hand, the Hamilton- Jacobi equation of (fl7|) 
can be rewritten as 



d t S- 



+ X^ + Xd q -v + q -(d q S + \^)-L = Q. (42) 



Applying Eq. (1411) . the above equation thus becomes 



d t S + H(<l,P)\p=d q s = 0, 



(43) 



q(q,d q S;t) 



dH 
dp 



p=d q S 



(46) 



Inserting Eq. (|4"6"|) into the constraint of Eq. (|2"TT) . one 

thus has 



d t n + d q ■ (p?= 



p=a q s 



0. 



(47) 



It is then imperative to check whether the above equa- 
tion is consistent with Eq. (j45]). To see this, notice that 
v in Eq. (|46|) is obtained by averaging q(±X). Hence, 
Eq. (|47l) has to be obtained by averaging Eq. (|4"5j) for 
the case of ±A as well. One can evidently see by tak- 
ing into account Eq. (fl9|) that for classical Hamiltonian 
at-most-quadratic in momentum, this is indeed the case. 
Namely, the constraint of Eq. (|4T|) is indeed consistent 
with Eq. ([4"5]). 

Given the classical Hamiltonian, we have thus two rules 
of Eqs. (|4l"|) and (|46|) to be applied to the classical me- 
chanical equations of (1431) and (|4"5j) and proceed in the 
way described in the previous section to arrive at the 
Schrodinger equation with a unique Hermitian quantum 
Hamiltonian. The above derived rules are just the rules 
of quantization proposed heuristically in Ref. [lj, that 
is Eq. (5) of Ref. [l|, where formal "replacement" there 
is re-interpreted in the present paper as physical "sub- 
stitution" [ll|. Hence, we have given a justification of 
the rules postulated in Ref. [l[ in term Hamilton- Jacobi 
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condition with a random constraint. In other words, the 
hidden random variable A postulated in Ref. [l[ is given 
physical interpretation as a random Lagrange multiplier 
that arises in Hamilton- Jacobi theory with a randomly 
fluctuating constraint. 

V. CONCLUSION AND DISCUSSION 

We have proposed a statistical model of quantization 
given a classical Lagrangian by assuming the existence of 
hidden random variables and accordingly imposing the 
Hamilton- Jacobi condition for the actual trajectory to 
a random constraint due to the fluctuations of hidden 
variables. Quantum fluctuations is shown to be emergent 
corresponding to a specific constraint depending uniquely 
on the choice of the classical Lagrangian and assuming 
that the Lagrange multiplier that arises in the dynam- 
ical system with constraint, which is fluctuating due to 
the fluctuations of the constraint, can only take binary 
values ±H/2 with equal probability. Given a classical sys- 
tem, the model leads to a unique quantum system with 
a straightforward physical interpretation. 

Let us mention several interesting problems that can 



be raised within the present statistical model. First, it is 
imperative to ask if the model can suggest new testable 
predictions beyond quantum mechanics. Such prediction, 
obtained by allowing fluctuations of |A| around H/2 with 
very small yet finite width, is reported in Ref. |l4|. Next, 
while we have offered an explanation on the physical ori- 
gin of Planck constant in term of Lagrange multiplier, 
there is still a missing explanation on what determines 
its numerical value. Such a question of course beyond 
the standard quantum mechanics and is permissible only 
within a model in which quantum fluctuations is emer- 
gent [lij]. Finally, notice that in the model, the ran- 
dom velocity constraint v is uniquely determined by the 
choice of classical (unconstrained) Lagrangian. In this 
sense, the random constraint is already inherent in or 
self-generated by the system being constrained. Hence, 
to each quantum system there is a hidden context specific 
to the system. 
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